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Abstract 

We investigate the generation of gravitational waves due to the gravitational insta- 
bility of primordial density perturbations in an early matter-dominated era which could 
be detectable by experiments such as LIGO and LISA. We use relativistic perturba- 
tion theory to give analytic estimates of the tensor perturbations generated at second 
order by linear density perturbations. We find that large enhancement factors with 
respect to the naive second-order estimate are possible due to the growth of density 
perturbations on sub-Hubble scales. However very large enhancement factors coincide 
with a breakdown of linear theory for density perturbations on small scales. To pro- 
duce a primordial gravitational wave background that would be detectable with LIGO 
or LISA from density perturbations in the linear regime requires primordial comov- 
ing curvature perturbations on small scales of order 0.02 for Advanced LIGO or 0.005 
for LISA, otherwise numerical calculations of the non-linear evolution on sub-Hubble 
scales are required. 

1 Introduction 

Gravitational waves are a probe of the very early universe that go beyond electromagnetic 
signals such as the cosmic microwave background. Gravitons can propagate essentially un- 
scattered from any energy scale below the Planck density. As a result they have been 
considered as probes of violent events in the early universe, such as bubble collisions [1], 
preheating after inflation [2, 3, 4, 5, 6, 7, 8, 9] or cosmic strings [10]. 

Much of the work to date has been based on the local generation of gravitational waves, 
e.g., due to time-varying quadrupole moment in flat spacetime [11]. But recent [12, 13, 
14, 15, 16, 17, 18, 19] work has developed cosmological perturbation theory to deal with 
the generation of gravitational waves at second order from first-order density perturbations. 
Such a relativistic treatment is required for inhomogeneities on or above the Hubble scale in 
an expanding cosmology. It is possible to make a perturbative calculation of the gravitational 
waves inevitably generated at second order from the existence of linear density perturbations 
on all scales. 



In the early radiation-dominated era, Ananda et al [17] found that gravitational waves 
are produced when density perturbations, which are overdamped at early times on super- 
Hubble scales {k/a <^ H) come inside the Hubble scale {k/a = H) and begin under-damped 
oscillations. The Hubble damping leads to the rapid decay of the amplitude of sub-Hubble 
metric perturbations, and an almost scale-invariant spectrum of gravitational waves is left 
on sub-Hubble scales to propagate freely, redshifted by the cosmological expansion. 

Baumann et al [18] extended this calculation to follow the evolution of the gravitational 
waves generated at second order through into the matter-dominated era, and ultimately 
the present late-time acceleration using the numerical solution for the evolution of linear 
density perturbations. Their work showed a surprising behavior in the matter-dominated 
era, where tensor metric perturbations grow on large-scales until reaching a constant value 
once they come inside the Hubble scale, producing a larger amplitude on scales close to 
the Hubble scale today than the first-order gravitational waves generated by inflation [18]. 
Although intriguing, it seems very difficult to detect such extremely long-wavelength tensor 
perturbations which are only produced at late cosmic times [13, 20]. 

In this paper we will investigate the production of gravitational waves in an early matter- 
dominated era, preceding the standard radiation-dominated era, before primordial nucleosyn- 
thesis. A matter-dominated era is expected to occur in the very early universe after a period 
of inflation driven by overdamped scalar flelds. At the end of inflation the Hubble damping, 
if, decreases and the scalar flelds become massive (m > H). Oscillating massive scalar flelds 
have an effectively pressureless equation of state [21]. The decay of the oscillating flelds 
leads to the reheating of the universe and the start of the standard Hot Big Bang model 
[22]. If the decay is slow F <S if at the end of inflation then we have an extended early 
matter-dominated era. An early matter-dominated era might also occur if weakly-coupled 
massive scalar flelds (or moduli) with non-zero vacuum expectation values come to dominate 
over the inflaton decay products sometime after inflation. In the curvaton scenario it is a 
weakly-coupled massive fleld, rather than the inflaton, whose inhomogeneous perturbations 
give rise to the primordial density perturbation [23, 24, 25, 26]. Gravitational waves gener- 
ated at second-order have previously been studied in the curvaton scenario [27], in the limit 
where the curvaton fleld does not dominate the energy density of the universe, i.e., with no 
early matter era. 

There are few constraints on an early matter era since the comoving Hubble scale must be 
less that a few parsecs. The primordial density perturbations on such small scales have long 
since been erased by Silk damping, and a stochastic background of gravitational waves might 
be one of the few ways we can probe the primordial power spectrum on such small scales 
(though primordial black holes might be another [28]). Gravitational waves produced in an 
early matter era would have wavelengths much smaller than the Hubble size at primordial 
nucleosynthesis and thus, depending on the reheating temperature at the end of the early 
matter era, could be directly directed by gravitational wave detectors currently in operation 
or being planned. 

This paper is organised as follows. In section 2 we summarise the results we will need 
for the evolution of linear density perturbations in a matter-dominated era. In section 3 
we give the evolution equation for second-order gravitational waves in a matter dominated 
era and present the resulting power spectrum for tensor metric perturbations at the end of 
the matter era. This leads to an effective density of gravitational waves at the present day. 
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^GW,o, presented in section 4. In section 5 we discuss the detectability of the gravitational 
wave background in possible early matter eras, and we conclude in section 6. 



2 Density perturbations in an early matter era 

In this paper we will work in the longitudinal [29] (or Poisson [12, 19]) gauge where the 
perturbed metric is 

ds'^ = a2(^)[-(l + 2^)(Pri + [(1 - 2^)5^^ + 2F(^,j) + hij]dx'dx^] (1) 

where rj is conformal time, and the conformal Hubble rate is 7i = a' / a where primes denote 
a derivative with respect to conformal time. Fj describes transverse vector perturbations 
and hij describes tensor metric perturbations which are transverse and trace-free. 

Vector perturbations, Fj, must vanish at first order in the absence of any vector part of 
the fluid 3-vclocity, e.g., during a period of inflation dominated by a scalar field. First-order 
tensor perturbations propagate as a free field, but arc rapidly rcdshiftcd to large scales by 
an infiationary expansion. Gravitational waves on smaller scales are produced from initial 
quantum vacuum fluctuations, but their amplitude depends on the energy scale of inflation. 
In the following we assume any flrst-order vector or tensor metric perturbations are negligible 
and focus on the second-order tensor perturbations that are produced from flrst-order scalar 
perturbations which are known to exist on observable scales. 

The scalar metric perturbations $ and ^' are supported by density and (curl-free) velocity 
perturbations of a fluid. In the absence of anisotropic stress we require $ = ^ [29] , and the 
linear evolution equation for $ is given by [19] 

$"H-3W + (27i' + 7i2)$ = 47rGa^5P (2) 

In a matter-dominated era with negligible pressure, SP — 0, we have a oarf and Ti. — 2/77, 
so the linear evolution equation for $ is simply 

$" + -$' = 0, (3) 

with the general solution ^ — C + D/rf , where C and D may be spatially inhomogeneous, 
but are constant in time. Considering only regular solutions at early times (77 — > 0) requires 
D = and we then have $ constant in time on all scales. 
We will often flnd it convenient to use a Fourier transform 

$(x) = -^ / rf^k^ke^''-, (4) 

where the initial power spectrum of an isotropic distribution of scalar metric perturbations 
will be a function of A; = |k| 

(^k$k') = ^'^'(k + k')TO (5) 
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Note that the hnear comoving density perturbation grows in time on all scales during a 
matter-dominated era and is given by [30, 19] 

^ = A.V^$. (6) 

where = S'^Widj is the comoving spatial Laplacian. Because we have assumed pressure 
is negligible, gravitational instability leads to a growing density perturbation on all scales. 
Eventually this will lead to a breakdown of the linear evolution on sub-Hubble scales when 
the density perturbations become of order one, corresponding to 

kNiiv) ~ v-^/^n^ > . (7) 

The above is an idealisation of any realistic model. In particular an oscillating mas- 
sive scalar field with mass m > H has a Compton wavelength A ~ and hence can 
only be described as pressureless matter on comoving scales k ^ am. Fourier modes with 
k > am correspond to relativistic modes with non-negligible pressure and we expect $ to 
be suppressed on these scales. For simplicity, and to remove any dependence on the pre- 
ceding cosmological evolution, we will assume in our calculations that there are no density 
perturbations on sub-Hubble scales at the start of the early matter era 

V{k) = for A; > A;dom (8) 

where k^om = ^dom < om is the comoving Hubble scale at the start of the matter-dominated 
era^. Any density perturbations on smaller scales will provide an additional, source for 

gravitational waves. 

Even so, density perturbations on our smallest scale will become nonlinear, fcdom > 
kNiirj)-! if the matter era lasts long enough. In what follows we will present results both 
extrapolating the linear results into the nonlinear regime, and also from imposing an abrupt 
cut-off on the scalar power spectrum at the non-linear scale, such that 

V{k) = for > fccut (9) 

where fccut = iiiin[A;dom, ^ArL(^)], which we expect to provide a conservative lower bound on 
the amplitude of gravitational waves generated at second-order. 

On very large scales the power spectrum of the primordial scalar perturbation is com- 
monly approximated by a power law 

where the numerical factor 9/25 comes from the relation between scalar curvature perturba- 
tion in the longitudinal and comoving gauges on large scales in a matter dominated era [29]. 
Observations of the cosmic microwave background (CMB) radiation [31] give the amplitude, 
K, 2.4 X 10~^, and spectral index, ^ 0.96, for the primordial density perturbations 
on very large scales 27rao//i;* ~ 100 Mpc today. We should note however that the primordial 
density perturbations could be very different on the much smaller scales relevant for the 
direct detection of gravitational waves, e.g., 2Trao/k ~ 1 km. 



^In the late universe this corresponds to the turnover in the matter power spectrum on scales k > fceq 
where fceq = Heq is the Hubble scale at matter-radiation equality. 
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3 Generation of tensor perturbations 

The existence of first-order scalar perturbations inevitably leads to second-order vector and 
tensor perturbations [13, 15]. 

Analogous to Eq. (4) for scalar fields, we can give the Fourier transform of the tensor 
metric perturbations ([17, 18]) 

-— Te'''"[^k(77)ei,(k) + h^eijik)] , (11) 

where e'^^{k) is the tensor polarization. The two polarization tensors eij{k) and eij{k) can be 
calculated in terms of the orthonormal basis: 

eij(fc) = -^[ei{k)ej{k) - ei{k)ej{k)] 

e,,{k) = -^[e,{k)e,{k)+e,{k)e,{k)] (12) 

where e and e are orthonormal transverse vectors, e.k = e.k = e.e = and e.e = e.e = 1. 
The power spectrum of the tensor perturbations is given by 

1 27r2 

{hMhM) = -^S'{k+k')rH{k,v), (13) 

where the factor of 1/2 arises since, by convention, Vh includes the contributions from both 
polarisations (/ik as well as hk). 

The evolution of second-order tensor mode in (1) is given by the wave equation [17]: 

/i^;, + 2m^ + = , (14) 

where the Sjj'" is a transverse and trace-free source term. If we include terms up to second 
order in the scalar perturbations, then Sj^'^ is the transverse and tracefree part of [17, 18] 



1 



3{l + w)W 

2c2 



ai(^' + 7i$)6>j(*' + 7i$) 



[37i(7i<l>-^') + V'^]ai6>j($-*) (15) 



where w — P/p is the equation of state and = P'/ p' is the adiabatic sound speed. 

In the approximation of a pressureless matter-dominated era we have -u; = = and 
$ = = C(x) and the source term simplifies considerably. If we substitute Eq. (4) into 
Eq. (15) we obtain 

S.j{x) = J d'^kd^'k' [-Ak'ik', - '^-k'~k,] e^(^+'^')-^ (16) 
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Substituting Eqs. (16) and (11) in Eq. (14) we find the evolution of the amplitude of each 
tensor mode in Fourier space during the matter era 

K + + k'^K = , (17) 



where 



= f (27r)-i j d3ke(k,k)*^_i^*i^, (18) 

and _ 

e(k, k) = e'^{\<i)kikj = sin^ 9 , (19) 

where 9 is the angle between k and k, 

cose = ^. (20 
kk 

3.1 Evolution of the tensor mode 

A striking feature of a matter dominated era is that the source term 5'k in the wave equation 
(17) is constant for linear density perturbations. This contrasts with, for example, a radiation 
dominated era, where and hence the source for tensor modes, decays on sub-Hubble scales 
[17]. 

The general solution of the evolution equation (17) in the matter era is thus 



A;2 Y 2k^r]^ 



Ci,J3{kr]) + D^Ysikri) 



S'k , „ f sm{ki]) - kr) cos{kr])\ f cos{kr]) + krj sin{kr])\ 

- F + ^'^'v J" 'V ) ^ ' 

where J and Y are the Bessel functions of the first and second kind and Ck and are 

constants of integration. 

The J and Y modes describe the damped, but source-free oscillations of the gravitational 
field and thus have the same form as the usual solution for first-order gravitational waves. 
The Y mode is singular as — * and rapidly decays at late times. J is regular at early 
times and gives oscillations whose amplitude redshifts with the expansion of the universe, 
|/ik| oc a~^. However the constant second-order source term, S, supports a constant tensor 
part of the metric perturbation at late times. This behaviour is quite different from that 
usually associated with gravitational waves, reflecting the fact that this is no longer a freely 
propagating wave, but rather a metric distortion sourced by terms quadratic in first-order 
scalar perturbations. 

If we impose the initial condition h — h' — when 77 = 0, then Ck = — SSk/Zc^ and the 
singular term is absent, Dk = 0, so we have the particular solution 



^^^^ 



1 + 3 



krj cos{kr]) — sm{kr]) \ 
k^rf J 



(22) 
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At early times, or equivalently in the large scale, super- Hubble limit (/c -C 7Y), we find a 
growing tensor perturbation 

= ^V' . (23) 

We see that the large-scale tensor mode grows at the same rate as the comoving density 
perturbation (6) in the matter era 

hoc-^cc-. (24) 

The tensor amplitude grows until the mode enters the Hubble scale and at late times, or on 
sub- Hubble scales {k ^ Ti), it becomes constant 

h^ = ^ (25) 

Substituting Eq. (18) for into Eq. (22) we will write the solution for the tensor per- 
turbation as 

where the growth function for the tensor modes is given by 



g(krj)^l + 3(^^^l^^^^^^ , (27) 

which approaches unity at late times on sub- Hubble scales. 

3.2 The power spectrum of the gravitational waves 

Prom (26) we can immediately write down the two-point function for the tensor modes 

(hMhr^^iv)) = (||[^)'^"'/'"'" d'ke{k,k)d'k'e{k',k') (28) 

where on the right-hand-side for a Gaussian distribution of scalar perturbations we have (for 
non-zero k and k') 

i^^-i ^k'-k' ^k') = (^k-E ^k'-E') (^k ^k') + (^k-k ^k') (^k ^k'-k') (29) 



Substituting Eqs. (29) and (5) into Eq. (28) we obtain 



V 3 y 2k^ 
X jd^ke{k,k) e(k',k) +e(k',k-k) :P(k - k) :P (k) , (30) 



where e(p, q) = g^sin 6 is defined in Eq. (19). 
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For simplicity we will assume that the power spectrum of the primordial scalar pertur- 
bations given in Eq. (10) is effectively scale invariant = 1 on scales k < /cdom- The power 
spectrum (13) of gravitational waves generated is then 

V,{k,n) = 2 (^^MMy (^^) h{k/kaom) , (31) 

where the integral 

h{k/kdom) = Tr4— I c^'k ~}^'^^~\ 0{kdom - k) e(A;dom - |k - k|) , (32) 



can be written as 



1 /•! (^ ,,2^2^,3 



Iiix)= I dn I dy — jQ{l- -y' + 2xyfi) , (33) 

J -I Jo (x^ + — '^xyjj) 2 

and we introduce the Heaviside step function, ©(/cdom ~ ?); to cut-off the scalar spectrum, 
P(g), on small scales q > /cdom- 

In the limit x <^ 1 the step function is equal to one throughout the integral in Eq. (33) 
and we obtain the analytic expression 

w N 16 4 16 2 

hix) ^ x^ . (34) 

^ ^ 15 3 35 ^ ^ 

For any x < 1 we can take Ji to be a numerical factor < 16/15. 

Equation (31) thus gives a simple numerical estimate of the tensor perturbation generated 
at second order during a matter dominated era. In particular at the end of an early matter- 
dominated era we have 



24\' . . /k. 



Vh{k, r^dec) ^ 2 [^-j X X g^k/k^,,) x h{k/kdor.) , (35) 

where k^ec = ''^dec the Hubble scale at the start of the radiation-dominated era, when the 
matter decays into radiation. This power spectrum for tensor perturbations is shown by the 
solid line in Figure 1 for an example where kdec = 10^A;doin- 

In the super-Hubble limit the tensor amplitude grows and we have g{kri) ^ k^^rf /IQ from 
Eq. (27). Thus on scales larger than the Hubble size at the end of the matter era {k < k^ec) 
we have 

n(A;,r^dec)-0.5A^(^) , (36) 

where we have taken Ii{k/kdom) — 16/15 using Eq. (34). The tensor perturbations thus have 
a steep blue spectrum on large scales, and are strongly suppressed on super-Hubble scales 
at the end of the matter era. 

On scales which enter the Hubble scale during the early matter dominated era k^ec < 
k < kdom then we find that the tensor amplitude has settled down to a constant value by 
the end of the matter era, g{kr]dec) — 1 and we have 

Vnik, 77dec) ^ 46A^ (^^^ /i(fc/fcdom) , (37) 
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Figure 1: The power spectrum of gravitational waves, shown as a function of wavcnumber k, 
generated from scalar perturbations during a matter dominated era, Vh{k,ri^cc)- The solid 
line shows the prediction using the linear matter power spectrum down to kdom = lO'^/cdeo 
the comoving Hubble scale at the start of matter domination. The dotted hne shows the 
prediction when the matter power spectrum is truncated at /Ccut — 200/cdec- ^dec denotes the 
Hubble scale at the end of the matter era. 
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Thus we find a decreasing power spectrum for tlie tensor modes on small scales with Vh oc k 
if we take Ii{k/kdom} to be constant. In practice Ii becomes small for k ~ k^om leading to 
an additional suppression on the smallest scales. 

Thus wc find the maximum amplitude of tensor perturbations is generated at the scale 
just entering the Hubble scale at the end of the matter era, kdec, for which we have 

V^^ = n(fcdec, %ec) ^ 50A^ (^^^ , (38) 

where we have taken /i(A;dec/^dom) — 16/15 for k^ec ^ ^dom using Eq. (34). 

Although the tensor power spectrum generated at second order is necessarily proportional 
to the square of the first-order scalar power spectrum, we find that the constant source 
term from scalar perturbations, Sij in Eq. (16), extending to sub-Hubble scales, leads to a 
additional factor k^om/kdec which may be large depending upon the duration of the early 
matter dominated era. Thus we find that the tensor power spectrum may be significantly 
enhanced with respect to the naive estimate Vh ~ A^. 



3.3 Nonlinear cut-off 

The gravitational wave power spectrum (38) becomes largest when the comoving Hubble 
scale at the end of the matter era, fcdcc, becomes much larger than the smallest scale fcdom- 
But as the Hubble scale grows, the comoving density contrast (6) becomes large on scales far 
inside the Hubble scale, signalling a breakdown of the finear results used thus far to estimate 
the source term in Eq. (17). 

Below the non-linear scale, k^L defined in Eq. (7), a pcrturbative analysis suggests that 
power will be rapidly transfered to smaller scales [32] leading to a suppression of $k and 
thus the source, 5*15;, on scales k > k^L- A full analysis of the nonlinear regime requires a 
full numerical treatment, such as a lattice field theory calculation, as has been performed 
in preheating models at the end of infiation, or an N-body simulation in the (small-scale) 
Newtonian regime. 

In practice we will obtain a conservative lower limit on the amplitude of primordial 
gravitational waves by assuming the scalar power spectrum is vanishing on all scales k > 
kNiiv) instead of the fixed comoving cut-off, fcdom- This implies a time-dependent source 
term (18) since the upper limit of the integral in k-space becomes time-dependent. 

For A; < /c^^ we find S^k^iV and thus 

V = W^^k-H^k. (39) 

The rate of change of the source term is thus slow compared with the decay time for the 
transient part of the solution in Eq. (21) for k > Ti and thus we take the quasi-static 
generalisation of Eq. (25) for sub-Hubble scales 

ft^M^M;^. (40) 

Thus a very conservative lower bound on the power spectrum of tensor perturbations 
generated on sub-Hubble scales at the end of an early matter-dominated era {kdec < k < /Ccut) 
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is given by the generalisation of (37) 



n(A;,r7dec)-46A^ 



(41) 



where we have k, 



;cut = min{kdoni, kMLiVdec)} and 



'dec 



200 A;, 



'dec ■ 



(42) 



For k ^ kNL we find that either or $k-k integrand in Eq. (18) vanishes for 

any k and the source term, S^, goes to zero. Equation (21) for ~ reduces to the 
standard solution for a free gravitational wave in a matter-dominated era, whose amplitude 
is redshifted, \hk\ oc a~^, on small scales, k > k^L ^ Ti-. Thus we obtain 



Hence Vh{k,ridec) oc /c"'', and the power spectrum of the gravitational waves is strongly 
suppressed on small scales {k ^ kNiirj)) if non-linear evolution suppresses the scalar per- 
turbation on these scales. 

Assuming the scalar perturbations rapidly decay to effectively zero on non-linear scales 
may be unduly pessimistic as in fact gravitational instability continues on small scales and 
non-linearity transfers power to smaller scales. However eventually the scalar metric per- 
turbation must decay on the smallest scales, where the velocity of the matter becomes non- 
negligible. We leave a detailed numerical calculation for future work and henceforth present 
the predicted gravitational wave background using both the linear result on all scales and 
the linear result cut-off at the non-linear scale. 

4 Present density of gravitational waves 

The effective energy density of a stochastic background of gravitational waves, on scales 
much smaller than the Hubble scale, is given by [33] 



Note that our second-order tensor modes (22) produced by linear scalar perturbations become 
constant on sub- Hubble scales in the early matter dominated era and therefore have negligible 
energy density at that time. They behaves quite differently from the conventional view of 
gravitational waves. But as the Newtonian potential $ decays on sub-Hubble scales either 
due to non-linear evolution on small scales or in the radiation era on sub-Hubble scales, the 
source term for tensor modes also decays on sub-Hubble scales and we are left with freely 
oscillating gravitational waves. 

The fraction of the critical energy density in gravitational waves per logarithmic range 
of wavenumber k in the radiation era is 




(43) 




(44) 




(45) 
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Figure 2: The present energy density of gravitational waves, ^gw,q, generated during a mat- 
ter dominated era shown as a function of wavenumber k. In this example F — (fcdom/ ^dec)^ — 
10^. The solid line shows the result predicted using the linear matter power perturbation for 
k < kdora, while the dotted line shows the result using the matter power spectrum truncated 
at k > kcut = SOO/cdec- 

During and after the radiation-dominated era, the density of gravitational waves on sub- 
Hubblc scales then redshifts exactly as any non-interacting relativistic particles and in the 
present day we have 



where the present density of photons is fl^^ ~ 1.2 x 10^^, and we neglect additional numerical 
factors due to the detailed thermal history, such as the heating of photons by the annihilation 
of other relativistic particle species. 

4.1 Linear scalar perturbations 

If we take Eq. (37) for the amplitude of tensor perturbations for k < kdec at the start of the 
radiation era, when Ti — k^ec, we have 



and this remains constant (assuming no further production of gravitational waves on sub- 
Hubble scales) during the radiation era. The present day density of second-order gravita- 
tional waves produced due to first-order scalar perturbations is thus given by 




(46) 




(47) 




(48) 
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The density as a function of wavenumber, k, is shown in Figure 2. 

Whereas the power spectrum (37) at the end of the early matter dominated era has a 
maximum on the Hubble scale at the start of the radiation era, /cdcc, we find that the present 
density of gravitational waves is largest on comoving scales of order the Hubble size at the 
start of the early matter era, kdom, 

/I \ 2 
A r\ I "'dom 



^Gw,o{f^dom) ~ ^n^y,o j ■ (49) 

We find that the maximum density of gravitational waves generated from linear density 
perturbations during an early matter dominated era is enhanced with respect to the naive 
expectation, ^GW,o{kdom) ~ ^^^7,0 ~ 3 x 10~^^, by a factor 

/, \ 2 

-i2 / "^dom 



V r>^dec / 

This enhancement factor represents enhanced amplitude with respect to the gravitational 
waves that would be produced by scalar metric perturbations Vh{kdoun ^dom) ~ for 
modes which re-entered the Hubble scale at the start of the early matter era but then 
redshifted once inside the Hubble scale. Instead wc find that the tensor perturbations remain 
constant even on sub-Hubble scales while they are supported by constant linear metric 
perturbations during the matter era, leading to 'P/i(fcdom, '7dec) ~ ^n- The energy density of 
such modes during the radiation dominated era is proportional to their frequency, A;^, and 

we find ^GW,o{kdoTn) ~ F^A^Q^^o- 

The enhancement factor calculated by continuing to use the linear result for scalar metric 
perturbations, $k = 0, down to very small scales, /cdom ^ ^dec, could be very large indeed. 
It is determined by the duration of the early matter-dominated era, for which Ti. — aH ex 
t~'^/^ oc if^/^, and we have 

F^=(^\'\ (51) 

If the early matter-dominated era arises due to a massive scalar field displaced from the 
minimum of its potential in the very early universe, then we require i/dom < ni-, the mass of 
the field, and the matter era will end when the field decays, ifdec ~ T, and thus 

< (^)^" . (52) 

In addition we require that the matter fields decay before primordial nucleosynthesis and 

thus r > i^BBN ~ McV^/Mpi where Mpj = G~^. In practice there is a tighter bound on the 
Hubble rate at decay for particles with mass m > 10 TeV, if we require that their decay 
rate should not be suppressed by more that the Planck scale, and thus T > m^/Mp^. In the 
optimal case where the field is precisely 10 TeV we find the weakest bound 



1!^ ^1020. (53) 



m J 



The bound is more restrictive for lighter fields, for which i^dom must be smaller, or more 
massive fields, which must decay earlier, unless their decay rate is suppressed with respect 
to gravitational strength decay. Nonetheless large enhancement factors would be possible in 
early matter-dominated eras that are sufficiently long-lived. 
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4.2 Non-linear cut-ofF 

In practice we have seen that the vahdity of hnear results for the scalar perturbations is 
restricted to a limited range of scales due to non-linear effects on scales k > kNL where 
^ATL is given by Eq. (7). Thus the maximum value of for which we can reliably use the 
calculation based on linear scalar perturbations is when kdom — ^ArL(?7dec) ^-nd hence 



If we assume that non-linear effects give a rapid suppression of the scalar metric perturba- 
tions on scales k > kNL then this non-linear cut-off suppresses the amplitude of gravitational 
waves, especially below the cut-off scale. For kNL < ^dom we find, from Eqs. (41) and (46), 
that for kdec < k < kNL 



The maximum value of the energy density of gravitational waves at the present time is 
reached for k ~ kNL'- 



For primordial density perturbations with the same amplitude on small scales as seen on 
CMB scales today, ~ 2 x 10~^, this gives ^owfii^NL) ~ 10~^^. This result neglects any 

source for the gravitational waves coming from non-linear scales, k > kNL, and thus gives 
a very conservative lower bound on the density of gravitational waves produced during an 
early matter era. 

5 Detectability 
5.1 Frequency ranges 

The comoving Hubble scale during a radiation-dominated era is given hj k = aH where 




(54) 




(55) 




(56) 




(57) 



where g^, is the effective number of degrees of freedom at temperature T. 
At the present time this corresponds to a physical scale 




(58) 



and a frequency 




(59) 
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The standard radiation dominated hot big bang must be restored before the epoch of 
primordial nucleosynthesis, and thus the temperature after the decay of an early matter- 
dominated era, Tdec, must be above 1 MeV and thus 

i/dec > 10"" Hz . (60) 

Note that ratio between the present frequency of Hubblc-scalc modes at the start of the mat- 
ter era and Hubblc-scale modes at the end of the matter era corresponds to the enhancement 
factor in Eq. (50) 

F=^. (61) 

^'dec 

Thus a large enhancement factor, F, also implies a wide range of frequencies over which 
gravitational waves are generated during an early matter era. 

Successful, standard, big-bang nucleosynthesis (BBN) [34] places an upper bound on the 
density of primordial gravitational waves, 

nGW,o<0.m^,oc^lO-\ (62) 

on any scales smaller than the Hubble scale at that time, corresponding to u > vbbn — 
10~^^ Hz. Gravitational waves generated during a early matter era are only present on 
LIGO scales, /ligo ~ 100 Hz, if Tdec < 10^° GeV, and on LISA scales, z^lisa ~ 10"^ Hz, 
if Tdec < 10^ GeV. Current limits from the Laser Interferometer Gravitational-wave Ob- 
servatory (LIGO) bound i^GW,o < 6 x 10~^ in the frequency range i^ligo ~ 100 Hz with 
Advanced LIGO sensitive to 1^gw,o ~ 10~^ in the future [35]. LISA should be able to de- 
tect a primordial GW background ^gw,o ~ 10^^^ at frequencies z/lisa ~ 10^^ Hz [36], and 
future experiments such as Big Bang Observer (BBO) may be able to detect a primordial 
background flcwfi ~ 10~^^ at frequencies z^bbo ~ 1 Hz [37]. 



5.2 Reheating after GUT-scale inflation 

Most inflationary models of the early universe incorporate an early matter dominated period 
immediately after inflation when the energy density of the universe is dominated by oscillat- 
ing, massive scalar fields. If the fields decay perturbatively to radiation this corresponds to a 
period of reheating. Current upper bounds from the CMB on first-order gravitational waves 
produced from vacuum fluctuations during inflation [38] place a bound on the maximum 
value of the energy density during inflation, = Minf < 10^^ GeV, and thus the Hubble 
scale at the start of the early matter era. 

Models such as chaotic inflation driven by a massive scalar field come close to saturating 
this bound and thus correspond to GUT-scale infiation. On the other hand constraints from 
the thermal production of gravitinos in supersymmetric models suggest that the maximum 
reheat temperature, Tdec, should be less than about 10^ GeV [34]^. This implies, from 
Eq. (51), a lower bound on the enhancement factor, > 10^, suggesting that the power 
spectrum of second-order gravitational waves generated during reheating after GUT-scale 

^Note that in the case of reheating at the end of chaotic inflation, the mass of the inflaton is required to 
be m ~ 10^^ GeV, and delaying reheating so that T^^c < 10® GeV then actually violates the bound given in 
Eq. (53) for this mass. 
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inflation could be larger than that from first-order gravitational waves at frequencies u ~ 
i^dom, and possibly much larger. 

Taking Eq. (10) for linear density perturbations on all scales k < k^om, we have, from 
Eqs. (49) and (51), 

where in the second line we set ~ 2.4 x 10~^. In this case the nucleosynthesis limit on the 
density of primordial gravitational waves (62) places a lower limit on the minimum allowed 
reheating temperature, T^^c > 10^ GeV, assuming Minf ~ 10^^ GeV. This value of O^owfl is 
obtained at high frequencies 

^.o--^xlO^H.(^)"'(^)'''. (65) 

which is always above the range of detectors such as LIGO. 

At the lower end of the range of generated frequencies, we have from Eq. (59) that 
z/dec ~ 10(Tdec/10'^ GeV) Hz, which is always below the LIGO range if Tdec < 10^° GeV. 
Using the linear perturbations (10) for all k < /cdom, the density of gravitational waves at 
LIGO frequencies is then, from Eq. (46), 

a™-,„(^,co) 5 X io»..-v3^^a„„ {J^y^ (S) 

Taking A|j ~ 2 x 10"^, this would be within the sensitivity of Advanced LIGO if the reheat 
temperature takes a low value, Tdec < 10^ GeV: 

10^ j (lo^^j (imI)' ('^^ 

However Eqs. (63-67) assume large enhancement factors, extrapolating the linear power 
spectrum (10) well into the non-linear regime on small scales, /cdom ^ defined in Eq.(7). 
If we assume non-linear effects cut-off the matter power spectrum below /catl — 200/cdec for 
~ 10~^, then from Eq. (55) we have for i/dec < J^igo < ^^atl 

n.„K,oo) - 50s.-/eA..a„„(^)"(ii^) (68) 

The amplitude increases for lower Tdec < 10^ GeV but it reaches a maximum value, ^gw,q ~ 
IQ-i^ at 

« 2 X Hz (j^) . (70) 
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and this drops below the LIGO range of frequencies for Tdec < 5 x 10 GeV. Advanced LIGO 
would thus only be able to detect a primordial gravitational wave background generated by 
density perturbations in the linear regime during reheating if the density perturbation ~ 
10~^ and Tdcc ~ 10^° GcV. However future experiments such as BBO could detect primordial 
gravitational waves generated during reheating even for ~ 10~^ if T^ec ~ 10^ GeV. 

5.3 Reheating after intermediate-scale inflation 

In many models of inflation the energy scale of inflation is set by the intermediate scale = 
Minf ~ 10^° GcV [39]. Such models include hybrid inflation [40, 41, 42] in which case the 
inflaton mass and Hubble-scale at the end of inflation is of order TeV. Tachyonic instability 
leads to the end of inflation when the inflaton reaches a critical value and the vacuum 
energy density driving inflation is rapidly converted to oscillating scalar fields. The spatially 
coherent oscillations of the inflaton fields rapidly fragment [43], but if a significant fraction of 
the energy density goes into weakly coupled massive fields which decay perturbatively then 
we may have an extended period of reheating after inflation. 

The lower energy scale during inflation, compared with GUT-scale inflation, leads to a 
lower frequency of gravitational waves produced at the end of inflation. Prom Eq. (65) we 
have 



If Tdcc < 10^ GeV then this is below the frequency range for LIGO but could be detectable 
by BBO if Tdec ~ 10^ GeV. Indeed if the inflaton or other moduli field excited at the 
end of inflation are weakly coupled and decay only just before big bang nucleosynthesis, 
^dec ~ 1 MeV, we have fdom ~ ^^lisa and a maximal enhancement factor ~ 10^^ for 
Minf ~ 10^° GeV, which would be large enough to violate the BBN bound in Eq. (62). 

Such large enhancement factors imply that the small scale density perturbations arc far 
into the non-linear regime. If non-linear evolution suppresses the scalar perturbations on 
scales k > kNL defined in Eq. (7), then the gravitational waves generated are suppressed for 
v > vnl where 



and we have set A^ 2.4 x 10 ^. Thus we find v^l < ^lisa for Tdcc < 10'^ GeV. For 
^dec ^ 10"^ GeV, LISA could detect this background if A^ ^ 10~^ on these scales. 

5.4 Moduli or curvaton domination 

Any scalar field, with mass m, displaced from the minimum of its effective potential after 
inflation begins to oscillate with an initial amplitude Xosc either at the end of inflation, when 
the Hubble rate is i^osc — -f^inf if m > i^inf, or once the Hubble rate drops below the mass. 
Hose — m, if m < ifinf. The fractional energy density of the field at this time is 




(71) 




(72) 




(73) 
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Note that ^x,osc < 1 since Xosc < Mpi , otherwise the field would slow-roll and drive a 
period of inflation. During a radiation dominated era following inflation, the density of the 
oscillating field grows relative to the radiation, eventually coming to dominate the energy 
density and drive an early matter era when the Hubble rate drops below [44] 

H,o^ ^ O^_i/osc ^(^X(^X(^Xrn. (74) 



x,osc-u». - y 3 y yH,,J \Mp,, 

The duration of the early matter era depends on the decay time of the field, i?dec = T, 
which must be before primordial nucleosynthesis. The duration is described by the factor 
defined in Eq. (51) 

F2 ^ 1020 -1/3 ( X^Y' ( J!lV ( ^ ( I}^y'^' (75) 

which can be very large if Xosc ~ -^Pi, becoming largest for fields which decay just before 
nucleosynthesis, T^cc 1 McV. Equation (75) is consistent with Eq. (63) for the amplitude 
of gravitational waves produced from reheating at the end of infiation if we set Xosc ~ -^Pi 
and Hose ^ m H^^f. 

If a moduli field is light during infiation, m < i?inf, then it acquires a spectrum of 
perturbations on large, super-Hubble scales, from initial vacuum fiuctuations on small, sub- 
Hubble scales. In the curvaton scenario [23, 24, 25, 26] these inhomogeneities about the 
background expectation value of the field is supposed to give rise to the primordial density 
perturbation, and hence in the simplest curvaton model in which the curvaton comes to 
dominate the total energy density of the universe before decaying, observations require [25] 



.37rXo; 

This determines the initial amplitude of the curvaton oscillations in terms of the infiationary 
energy scale, Xosc ~ A^^ifinf, and hence from Eqs. (49) and (75) we can calculate the peak 
amplitude of gravitational waves generated from linear density perturbations 

- ^ X (ra)"" iwkvf iM"^ • 

where we have taken ~ 2.4 x 10~^ and m = Hose < -f^inf , and we note that there is a lower 
bound on the Hubble rate during infiation, Hi^i > 10^ GeV [45] in the simplest curvaton 
models. Note also that if the curvaton comes to dominate the energy density of the universe 
before it decays, the primordial density perturbation is close to Gaussian (the non-linearity 
parameter /nl ^ -5/4 [46, 47, 48, 49, 50]). 

Thus we see that an early matter dominated era may have a long duration in the curvaton 
scenario leading to a large enhancement factor in our calculation based on the linear density 
perturbation. However, as in the case of reheating after intermediate-scale infiation, the 
frequency is below LIGO scales if the decay temperature is low. The highest frequency GW 
generated in the curvaton-dominated era would be at frequency 
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where the lowest frequencies would be, from Eq. (59), 



If non-linearities suppress the density perturbation on sub-Hubble scales in an early 
matter era then the only gravitational waves that are not reshifted away will be on the 
lowest frequencies, v < z/^l and we recover the same results as in the case of reheating after 
inflation, given in sub-sections (5.2) and (5.3) determined by the final decay temperature, 

T 
dec- 



6 Conclusions 

In a matter-dominated era the behaviour of the linearised density perturbations is particu- 
larly simple. The gravitational instability of pressureless matter leads to the growth of the 
comoving density perturbation, while the Newtonian metric potential, $, remains constant 
on all scales. This provides a constant source term, .Sk, at second-order in the wave equation 
for tensor metric perturbations, leading to the general solution given in Eq. (21). At late 
times this approaches the constant solution, given in Eq. (25), on sub-Hubble scales. 

This constant tensor mode during the matter era is quite different from the behaviour 
we would expect for a free gravitational wave on sub-Hubblc scales. Indeed the free part 
of the tensor perturbation, given by the oscillating terms in Eq. (21), are redshiftcd away 
|/ik| oc a~^, diluting any pre-existing gravitational waves on sub-Hubble scales at the start 
of a matter-era. But the constant amplitude of tensor metric perturbations in the Poisson 
gauge, supported by first-order metric perturbations in an early matter era, oc S^/k"^, 
becomes the initial condition for the amplitude of freely oscillating gravitational waves in 
the subsequent radiation era when $ rapidly decays on sub-Hubble scales. 

We have seen that a constant amplitude of scalar metric perturbations, ~ A^, ex- 
tending down to sub-Hubble scales can lead to a large enhancement in the amplitude of 
gravitational waves. On the smallest scales which we considered, fcdom corresponding to the 
comoving Hubble scale at the start of matter domination, the enhancement factor is easily 
understood due to the amplitude of tensor perturbations on this scale remaining constant 
during the matter dominated era instead of being redshifted on sub-Hubble scales. As a 
result the energy density of gravitational waves in the late universe is enhanced with respect 
to the naive estimate ^Igw,o ~ ^^^7,0) by a factor = (A;dom/^dcc)^i where /cdcc is the co- 
moving Hubble scale at the end of the matter era. For F"^ > rA^ ~ 4 x lO^r, where r <^ 1 
is the usual tensor-scalar ratio at first-order [38], the second-order gravitational waves gen- 
erated from linear perturbations during an early matter dominated era would dominate over 
those produced directly from vacuum fluctuations in the gravitational field during inflation. 
Large enhancement factors occur when there is a sufficiently long matter dominated era in 
the early universe, preceding the standard radiation dominated era required for successful 
primordial nucleosynthesis. This is certainly possible given our state of ignorance of the 
thermal history of the universe before nucleosynthesis. The spectrum of gravitational waves 
could provide one of the few constraints on the universe before nucleosynthesis. 
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However, a sufficiently long-lived matter era also predicts a breakdown of linear perturba- 
tion theory for density perturbations on sub-Hubble scales. Predictions based on linear the- 
ory can only be trusted for < (/cATL/^dec)^ ~ ^tz- ^ constant scalar metric perturbation 
in linear theory in the longitudinal gauge, $, implies a growing comoving density perturba- 
tion, Sp/p. If non-linear evolution of the matter perturbations suppresses the scalar source 
term on scales k > k]siL then the peak in the gravitational wave density is shifted from /cdom 
to kj^L and the maximum amplitude at this peak is reduced to ^gw,o ~ SA^Q^ q ~ 10"^^ 
for ~ 2 X 10~^ which would only be detectable by experiments in the far future such as 
the Big Bang Observer. 

In summary wc have seen that growing density perturbations during an early matter 
era could lead to a detectable background of primordial gravitational waves. But this is 
only detectable with planned experiments such as LIGO and LISA if the primordial metric 
perturbations are significantly larger on small scales (A^ ~ 6 x 10~^ for LIGO, or A^ ~ 
3 X 10~^ for LISA) than currently observed on very large scales, where A^ ~ 2 x 10~^, or 
if significant non-linear density perturbations survive on sub-Hubble scales during the early 
matter era. This latter case requires non-linear (and hence probably numerical) solutions for 
the density perturbations. Because these perturbations are on sub-Hubble scales, Newtonian 
theory may be sufficient on these scales. There has already been investigation of gravitational 
waves generated by nonlinear density perturbations on sub-Hubble scales during preheating 
after infiation [2, 3, 4, 5, 7, 8, 9], and our work suggests that similar numerical calculations 
will also be required to investigate whether a detectable level of gravitational waves could 
be produced from gravitational instability in an early matter era. 
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